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1, INTRODUCTION 

Metaheuristics are the best algorithms used to solve the NP-Hard optimization problems for 
which there 1s no exact effective known method. Sometimes, they can be the only possible methods to find a 
good solution to these hard problems. The most important advantage of using a metaheuristic is that it uses a 
high abstraction level. Therefore, it can be applied to wide different problems. New metaheuristics are 
proposed recently such as the hunting search [1], the chicken swarm optimization [2], [3] and the golden ball 
algorithm [4]. 

Hunting Search (HuS) is a developed metaheuristic for solving continuous optimization problems. It 
starts to be adapted to solve combinatorial optimization problems such as The Traveling Salesman Problem 
[5], [6] and the no-wait flow shop scheduling [7]. It 1s also proposed as a first adaptation for the Quadratic 
Assignment Problem [8]. HuS is a method inspired by group hunting of some animals such as wolves that 
organize their position to surround the prey. Each of them is relative to the position of others and especially 
in relation to the position of their leader. 

Quadratic Assignment Problem (QAP) [9] is the well-known discrete optimization problem of the 
category of the facilities location problems NP-hard. The problem is to assign a set of facilities to a set of 
locations in order to minimize the total cost of assignments. The QAP has several applications in 
combinatorial optimization problems such as Backboard Wiring [10] and scheduling [11]. In order to solve 
the QAP, several metaheuristics have been proposed [12]. The present paper proposes an Improved Hunting 
Search algorithm ({[HuS) adapted as a combinatorial optimization method to solve the QAP in a minimum 
CPU run time. 
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This paper is divided into six main sections. The first one is an introduction of the given method and 
the benchmark problem. The second section presents the metaheuristic HuS. The third section provides a 
detailed description of the QAP. The fourth section presents the proposed adaptation and improvement of 
HuS for the QAP. Numerical results obtained by the use of IHuS on the QAPLib instances [13] are 
demonstrated in the five section, and the last section concludes the whole work. 


2. HUNTING SEARCH ALGORITHM 
Hunting Search algorithm is an evolutionary method inspired by cooperative hunting of some 
carnivores that hunt bigger and faster preys than themselves with fewer energy. It is a population-based 
stochastic metaheuristic. The population is the hunting group that contains a set of solutions of the studied 
problem; each hunter represents a solution. The leader represents the best solution defined by an objective 
function. A hunter is characterized by its position that defines the distance between it and the other hunters. 
The hunting process in nature represents the search of optimum in the algorithm. The movements of 
the hunters to encircle the prey represents the operations of improvement of the initial solutions. Hunters do 
three main movements. Two movements of intensification operations. These are the movement toward the 
leader and the movement toward the other hunters. The third one is a diversification operation the 
reorganization of the hunters when they become very close to each other. HuS algorithm is presented in 
Figure 1: 
- HG (Hunting Group) is the population of hunters used for the search of the optimum. 
- NE (Number of Epochs) is the number of loops to make for the search. 
- IE (Iteration per Epoch) is the number of times per epoch where hunters make a move. 
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Figure 1. Flowchart of hunting search algorithm 


The main operations of HuS are: 
a) Initialize the Hunting Group: generate a set of random solutions for a problem; each hunter represents a 
solution. 
b) Moving Toward the Leader: hunters move toward the leader, which makes them stay closer to the leader 
until finding a better solution and becoming the new leader. 
c) Cooperation among Hunters: each hunter moves toward the other hunters or just makes a random move. 
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d) Reorganization of hunters: at some point, hunters become very close to each other, which means that they 
represent almost the same solution. In this case, they have to be regenerated, except the leader. 


3. QUADRATIC ASSIGNMENT PROBLEM 
3.1. Definition 

QAP is a combinatorial optimization problem of the class NP-Hard whose computational 
complexity increases exponentially by increasing the number of facilities. No exact method can solve the 
problem when the number of facilities 1s upper then twenty. 

Given a set of facilities to assign to a set of locations, there is a required flow between every two 
facilities and a required distance between every two locations. The problem is to find the best assignment of 
the facilities to the locations to have the minimum total cost of flows and distances. 

Figure 2 is an example of assignment of three facilities to three locations. D1 and D2 are the 
distances between the locations, and F1 and F2 are the flows between the facilities. 
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Figure 2. Assignment of three facilities 


The optimal solution is defined by an objective function from a discrete subset of the feasible 
solutions. 

Let E be the set of the feasible solutions; S is a subset of E; and g:S — R 1s the objective function. 
The problem is to find: 


min{g(s):s € S} (1) 


Where s is a solution from S. It is a vector of N, which contains indexes of facilities assigned respectively to 
locations 1, 2,...,n. n 1s the number of the locations. 
The objective function gives the cost of the assignment, defined as follows: 


g(s) = » > ftosen x dij (2) 
L=1 1 


j= 
Where (fj) € IR and (d;;) € R are two matrices. (f;;) is the square matrix that represents the required flow 
between facility 7 and j. (dj;) is the square matrix that represents the distance between location i and j. 


s € S, s(i) 1s the location to which facility 7 1s assigned and nis the dimension of the two matrices. 


4. IMPROVED HUNTING SEARCH ALGORITHM 

HuS as a metaheuristics uses intensification and diversification operations. The move towards the 
leader and the cooperation among hunters are the intensification operations; and the reorganization of hunters 
is the diversification operation. HuS uses also a population of individuals (the hunters) in the search space 
that undergoes mutation operations. It is an evolutionary algorithm. 
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4.1. Adaptation 
4.1.1 Initialize the hunting group 

To represent a solution of QAP in a program, an array of integer is used, where the array indexes 
refer to the locations respectively /, 2, ..., n (n is the number of locations) and the array contents refer to 
facilities assigned to each location. Figure 3 is an example of a QAP solution that represents an assignment of 
five facilities to five locations; the solution is represented by a hunter. 

In this first step, a set of random solutions called the Hunting Group (HG) is created. Each solution 
is represented by a hunter and the size of the set of solutions is called the Hunting Group Size (HGS). 
Figure 4 is an example of a HG of HGS=3. 
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Figure 3. An example of an assignment facilities Figure 4. An example of a HG initialization 
solution 





4.1.2 The movements of a hunter 

The only possible representation of the movement of a given hunter in the program is the 
permutation of two of its array contents. The algorithm is used to define the strategy and the sizes of 
permutations to do. Therefore, a movement of a hunter is equal to a permutation as it 1s shown in (3). 


1 movement = 1 permutation (3) 

The movement of a hunter H; toward another hunter Hj 1s also a permutation in a given array index 
k. One search for the location of the content array in Hj that is similar to the content of array Hj in the array 
index k and permute them in Hj as it shown in (4). Hj is the updated hunter and Hj is from which the update is 
inspired. 


1 movement(H;, H;, k) = 1 permutation(H,, H;, k) (4) 


Figure 5 is an example of movement of the hunter Hi toward the hunter H3 in the array index k=2. 
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Figure 5. An example of a movement toward a hunter 
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4.1.3 Moving towards the leader 

After the initialization of the HG starts the loop of the NE iterations. In this loop starts the second 
loop of IE iterations, where comes the move towards the leader (MTL). In this step, each hunter moves 
toward the leader by copying a part from the best solution; the size of the movement towards the leader 
(SML) which is the size of the copied part is calculated as follows: 


SML = rand X MML x D(H,,H;) (5) 


Where: 

- Rand is a uniform random number that varies between O and 1. 

- MML (Maximum Movement toward the Leader) is a number between O and | representing the 
maximum closer rate of a hunter to the leader. 

- Function D refers to the distance between two hunters. It is the number of the different 
assignments in the two represented solutions. HL is the leader and Hi is the hunter number 1. 

Each movement toward a leader (MTL) is characterized by a start array index that is chosen 
randomly and a size of movements toward the leader that is SML. 


1 MTL = SML x permutation (6) 


Figure 6 1s an example of a movement of the hunter number one towards the leader, the SML=2. 
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Figure 6. An example of a movement towards the leader 


4.1.4 Cooperation among hunters 

Always in the loop of IE iterations, after the move towards the leader comes the cooperation among 
hunters, where each hunter moves towards the other hunters by copying a part from their solution with the 
probability HGCR, or by changing their position relatively to themselves with the probability (1-HGCR). 
HGCR (Hunting Group Consideration Rate) is a number between O and 1. Figure 7 is an example of the 
movement of the hunter number | towards three different hunters (the case of the probability HGCR). And 
Figure 8 is an example of changing the position of a hunter relatively to itself (the case of the probability 1I- 
HGCR). 
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Figure 7. An example of a movement towards others hunters (the case of HGCR) 
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Figure 8. An example of changing the position (the case of 1-HGCR) 


4.1.5 Reorganization 

The intensification operations reduce the distances among the hunters. So, the solutions presented 
by these hunters become almost similar, which can lead to a blockage in a local optimum. Then, a 
diversification operation is needed: the reorganization hunters. Figure 9 is an example of reorganization of a 
population of hunters of HGS=3. 
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Figure 9. An example of reorganization hunters 


4.1.6 Strategies 
a) Backup system 

At the end of every intensification operation, a backup system is applied and only the good new 
solutions are taken. If the movement towards the leader or towards other hunters gives worse solution, one 
returns to the previous position. This strategy helps a rapid convergence towards local optimum. 


b) End criteria 
The end of the search process is determined by the number of epochs NE or by its half in the case of 
a stationary optimum. 


4.2. Improvements 
4.2.1 One by one backup system 

The first improvement is applied at the strategy of the backup system. Therefore, a backup is made 
for every change during the operation of intensification and not until the end. Figure 10 is an example of the 
improved backup system during the movement towards leader. 


4.2.2 Leader moving towards hunters 

In the initial HuS algorithm, the leader does not move along the search. The new proposed 
improvement is to move the leader towards the other hunters (LMTH) while using the backup system for 
preserving its quality of the best-found position. Figure 11 is an example of a leader moving towards three 
hunters using the improved backup system. 
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Figure 11. An example of a movement of the leader towards the other hunters 


4.2.3 Dynamic parameters 

The final proposed improvement is at the level of changing parameters dynamically during the 
search of optimum. This improvement helps better control the intensification and the diversification 
operations. Parameters are classify in three categories: parameters of static value, parameters of dynamic 
value and parameters of hybrid value. The parameters of static value are HGS, HGCR and NE; they have the 
same value from the beginning of the search until the end. The parameter of dynamic value is MML,; it 
changes its value during the search. And the parameter of hybrid value is EPS; it changes its value under 
certain conditions. 

The proposed definitions for the static parameters are as follows: 


HGS = N (7) 
NE = 100 (8) 


Where WN 1s the size of the evaluated QAPLib instance. 
The proposed definitions for the dynamic parameters are as follows: 


_ (EN -TN) 
RLW = Sn 5) 
MML =0.2+0.3 x RLW (10) 
EPS = min{g(sy,) — g(s,): 5s € 8S} (11) 
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oe Ee +0.01X EPS, RLW>0.5 (12) 
~ LEPS—0.01x EPS, RLW <0.3 

- he +5x(EN—TN), RLW>0O.5 (13) 
~ UE-5x(EN-TN), RLW<0.3 


Where: 
1. RLW (Rate of the distance between the Leader and the Worst hunter). 


2. IN (Trapped Number) is the number of epoch that the distance between the leader and the worst hunter is 
lower than EPS. 

3. EN (Epoch Number) is the number of epochs reached during the search. 

4. EPS (Epsilon) is the minimum distance between the leader and the worst hunter. 

5. sis a solution from the subset of solutions S. 

6. g(s_,) is the value of the objective function of the best solution represented by the leader during the 


search of the optimum. 
7. g(Sy) is the value of the objective function of the worst solution during the search of the optimum. 

According to (10) the value of MML is decreasing while the value of TN is increasing and vice 
versa. The minimum value of MMT is 0.2 and the maximum value is 0.5. 

The value of EPS in (11) is determined as the minimum distance during the search, then, in (12) is 
an improvement of EPS. 

According to (13) if the REW is superior to 0.5, which means that the distance between the leader 
and the worst hunter is 50% superior to EPS, the value of JE is increasing to intnsify the search. And if the 
RLW 1s inferior to 0.3, which means that the distance between the leader and the worst hunter is 30% infirior 
to EPS, the value of IE is decreasing to reduce the intensification during the search. The minimum value of 
IE is 30 and the maximum value is 100. 


5. EXPERIMENTAL RESULTS 
5.1. Experimental results of the proposed improvements 

This section presents the performance of [HuS algorithm on instances of QAPLib. The tests are 
performed on a computer processor Intel(R) Core(TM) 15-4300 CPU @ 1.9GHz @ 2.50 GHz and 4 GB of 
RAM. 20 times tested for each instance. The most important collected data from the obtained results over the 
20 runs are: 
Obest: The Best-Found Solutions (BFS). 
dave: The average of the BFS. 
PSD: Percentage of the Standard Deviation. 
Suc.: The percentage of Success in getting the Best-Known Solution (BKS). 
Oave: The average percentage of error in getting the BFS. 
Tavg: The average run time of the program in getting the BFS. 
Dbest: The best run time of the program in getting the BFS. 

PSD 1s calculated as follow: 


ee 


SD 
PSD = oe) 


avg 


2 
a= Lisi (BFS; ~ Savg) (15) 
\ 20 


Where 

e SD 1s the Standard Deviation. 

e BSF;is the best-found solution in the test number 1. 
Qavg 1S calculated as follow: 


Savg — BKS 
02> (Ce) x 100 (16) 


x 100 
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5.1.1 One by one backup system 

Table 1 shows the used parameter values for the tests of the new proposed backup system on the 
QAPLib instance Bur26a. Table 2 shows the obtained results compared to the old backup system (backup 
system for all the modified part in the solution). 


Table 1. Parameter Values 
HGS MML IE NE 
26 0.3 30 100 


Table 2. One By One Backup System Compared To Full Part Backup System 


Backup system Suc. (%) Bave (%) I best (sec) T'avg (sec) 
Full part 1 0.2049 2.145 2.911 
One by one 87 0.0107 0.901 24.634 


According to the table, the success of finding the optimum by the new proposed backup system has 
increased considerably compared to the one obtained by the old backup system, exactly 86 times better. The 
percentage of error is reduced more than 20 times, and the best time is reduced more than 2 times. Finally, 
the average time is increased more than 10 times because of the additional tests. 


5.1.2 Leader moving toward hunters 
Table 3 shows the obtained results for the tests of the new proposed operation: the move of the 
leader toward hunters (LMTH) on the QAPLib instance Bur26a. 


Table 3. Results For The Lmth Operation 
Proposed operation Suc. (%) Bave (%) I best (sec) Tavg (sec) 
LMTH 100 0 0.954 16.669 


According to the table, the success to find the optimum by the new proposed LMTH operation has 
reached 100%, and the average time is improved by 33%. 


5.1.3 Dynamic parameters 
5.1.3.1 Hunting Group Size 

Figure 12 shows the obtained results from the tests of the proposed values of HGS applied to the 
QAPLib instance Bur26a. The proposed values of HGS are: 


HGS = N (17) 
HGS = Nx2 (18) 
HGS = N+2 (19) 


Where N is the size of the evaluated QAPLib instance. N(Bur26a) = 26 

According to Figure 12, (18) is better than (19) in all the obtained results. (17) is better than (18) in 
all the obtained results except for Vavg by less than 10%. The most important result is Oavg. So, (17) is the 
better value for HGS. 


5.1.3.2 Maximum Movement toward the Leader 

Figure 13 shows the obtained results from the tests of the fixed maximum and minimum values of 
MML applied to the QAPLib instance Bur26a. According to Figure 13, the average error is equal to 0 and 
the average time 1s the least when the MML value is between 0.2 and 0.5. 
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Figure 12. The results obtained by changing the value of HGS 
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Figure 13. The results obtained by changing the maximum and minimum values of MML 


5.1.3.3 Epsilon and the iteration per epoch 
Figure 14 shows the obtained results from the tests of the RL'W values that increase or decrease the 
values of Eps and IE applied to the QAPLib instance Bur26a. 
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Figure 14. The results obtained by changing the Eps and the IE values 


According to Figure 14, the average error is equal to 0 and the average time is the least when the 
Eps and the IE values are improved when the value of RLW is inferior to 0.3 or superior to 0.5. Table 4 
shows the obtained results for the tests of the dynamic parameters compared to the static parameters on the 
QAPLib instance Bur26a. According to the table, the dynamic parameters have improved the best time and 
the average time to find the optimum. 
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Table 4. Results For The Dynamic Parameters 


Parameters Suc. (%) Bave (%) I best (sec) Tavg (sec) 
Static 100 0 0.455 17.720 
Dynamic 100 0 0.317 14.639 


5.1.4. [HuS applied to QAPLib instances 
Table 5 shows the obtained results by applying [HuS to 36 instances of QAPLib. Among 36 


instances of QAPLib evaluated in the Table 5, the proposed IHuS has solved 20 instances with no error, 13 
instances with error and has not solved 3 instances. 


Table 5. Numerical Results Obtained By Ihus Applied To Some Qap Instances Of Qaplib 


N° Instance BKS davg Bave (%) T'avg (sec) PSD Suc. (%) Ibest (sec) 

1 bur26a 5426670 5426670 0 14.640 0 100 0.317 

2 bur26b 3817852 3817852 0 13.547 0 100 2.824 

3 bur26c 5426795 5426795 0 3.218 0 100 0.688 

4 bur26d 3821225 3821225 0 7.951 0 100 1.579 

5 bur26e 5386879 5386879 0 4.221 0 100 0.653 

6 bur26f 3782044 3782044 0 3.159 0 100 0.297 

a bur26g 10117172 10117172 0 15.829 0 100 2.031 

8 bur26h 7098658 7098658 0 3.842 0 100 0.923 

9 tailOa 135028 135028 0 0.542 0 100 0.049 
10 tailOb 1183760 1183760 0 0.197 0 100 0.037 
11 tail2a 224416 224416 0 0.640 0 100 0.073 
12 tail2b 39464925 39464925 0 0.717 0 100 0.058 
13 tail5a 388214 388291 0.0198 4.855 0.08433 95 0.531 
14 tail5b 51765268 51765268 0 0.415 0 100 0.089 
15 tail7a 491812 493948.7 0.4340 26.937 0.41487 40 2.037 
16 tai20a 703482 709926.3 0.9160 23.782 0.38367 10 6.083 
17 tai20b 122455319 122510738.4 0.0453 9.483 0.13570 90 0.305 
18 tai25a 1167256 1187017.9 1.6930 56.834 0.38947 0 49.033 
19 tai25b 344355646 344367602.6 0.0035 18.168 0.01513 95 0.883 
20 tai30a 1818146 1851329.6 1.8251 112.737 0.37674 0 96.424 
21 tai30b 637117113 637271686 0.0243 91.640 0.03810 20 11.169 
De tai35a 2422002 2477518.6 22922 204.162 0.49207 0 173.891 
23 tai35b 283315445 283487663.15 0.0608 118.488 0.08008 60 27.158 
24 tai40b 637250948 63728605 1.4 0.0055 176.889 0.01111 70 29.569 
25 tai64c 1855928 1855928 0 107.5595 0 100 33.971 
26 lipa30a 13178 13178 0 17.403 0 100 1.091 
ZY lipa40a 31538 31777.75 0.7602 280.597 0.46111 20 181.749 
28 lipa70b 4603200 4603200 0 529.561 0 100 279.521 
29 escl6a 68 68 0 0.0615 0 100 0.011 
30 esc32a 130 132.1 1.6154 123.5548 1.12026 2 41.088 
31 esc64a 116 116 0 1.6071 0 100 0.306 
32 esc128 64 64 0 507.5493 0 100 23.036 
33 had20 6922 6922 0 0.9235 0 100 0.141 
34 kra30a 88900 89279.5 0.426 57.006 0.58714 65 4.080 
35 kra30b 91420 91590 0.186 120.770 0.10775 5 17.204 
36 chr25a 3796 4201 10.669 57.264 4.24970 5 44.783 


5.2. Comparison with other metaheuristic 

Figure 15 compares average percentage of error in getting the BKS of ten QAPLib instances 
obtained by the proposed I[HuS and the existing Genetic Algorithm (GA) [14]. The adaptation of GA for the 
QAP is the only contribution found in literature of an evolutionary algorithm that does not use the local 
search to improve the results. Therefore, it is the only available work in literature to compare results with the 
proposed work. 

According to Figure 15, it is clear that the average percentage of error in getting the BKS obtained 
by I[HuS is much better compared to that obtained by the GA. Knowing that the GA has been run on a better 
computer with Intel(R) Core(TM) 17-3770 CPU @ 3.40GHz and 8.00 GB RAM. 


Improved hunting search algorithm for the quadratic assignment problem (Amine Agharghor) 


154 O ISSN: 2502-4752 


s 6 

oD 

= 

S 4 
e 
0 


tai20a tai20b tai25a tai25b tai30a tai30b tai35a tai35b tai40a tai40b 


mHuS g@ GA 


Figure 15. The average percentage of error in getting the BKS obtained by I[HuS and GA for some instances 
of QAPLib 


6. CONCLUSION 

In this work, a new Improved Hunting Search algorithm is presented. It is the best evolutionary 
algorithm proposed to solve the quadratic assignment problem. QAPLib is the used library of the QAP 
instances to assess the performance of [HuS. To show the quality of [HuS, no local search is used. Three 
effective improvements are proposed: the one by one backup system, the leader moving toward hunters and 
the dynamic parameters. Very good results are obtained by applying [HuS to 36 instances of QAPLib 
instances. Finally, it 1s shown that the results obtained by [HuS are much better to that obtained by the 
existing evolutionary algorithm in literature. 
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